Infinite order excitonic Bloch equations for asymmetric nanostructures 
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We present a new exciton-based formalism for calculating the coherent response of asymmetric 
semiconductor multiple quantum well structures to ultra-short optical pulses valid to infinite order 
in the optical field and including the self-generated intraband fields. We use these equations to 
calculate and explain the oscillations with time delay of peaks in the spectrally-resolved degenerate 
four wave mixing signals from biased semiconductor superlattices, obtaining good agreement with 
experiment. 
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Over the last decade, there has been a great deal of in- 
terest in the nonlinear optical response of semiconductor 
nanostructures to ultra-short optical pulses. Many phe- 
nomena have been treated using the dynamics controlled 
truncation (DCT) theory of Axt and Stahl 0, whereby 
the response is expanded in powers of the exciting laser 
field, and is formulated to include n-point correlations 
111 L4 LI ■ This theory has been very successful, and cor- 
rectly deals with the intra- and inter- excitonic correla- 
tions to arbitrary order in the optical field. However, it is 
intrinsically a perturbative technique, while the theoret- 
ical description of many interesting nonlinear processes 
requires a non-perturbative approach. 

The semiconductor Bloch equations (SBE's) in the 
Hartree-Fock (HF) approximation provide a very com- 
pact description of dynamics to infinite order in the opti- 
cal field 4] . They have been used to describe Rabi split- 
ting near zero detuning || and the ac Stark effect . It 
has been shown SSl> however, that the HF factoriza- 
tion of four-particle correlation functions into products 
of two-particle interband polarizations and electron and 
hole populations incorrectly predicts that intraband cor- 
relations lose any excitonic signature as time progresses, 
and act simply as free electron-hole pairs at long times. 

In recent years, there has been a great deal of work on 
the nonlinear optical properties of asymmetric coupled 
quantum well systems such as biased semiconductor su- 
perlattices (BSSL's) US SHE El El El- Of partic- 
ular interest has been the interplay of the interband and 
intraband response in two-beam degenerate four-wave 
mixing (FWM) experiments, where ultra-short optical 
pulses with wave vectors ki and separated by a delay 
time t 2 i impinge on a sample S0SSE1- The inter- 
band polarization is responsible for pump-probe, four- 
wave mixing and higher order signals, while the intra- 
band polarization describes the coherent motion of the 
electrons within each band. In a BSSL, the electron 
energy levels form the so-called Wannier-Stark ladder 
(WSL): E n = Eq + nfkds, where lub = eErjcd/h is the 
Bloch frequency, d is the superlattice period, and Erjc is 
the applied along-axis DC electric field. Excitation of a 



coherent superposition of WSL states results in oscilla- 
tions in the intraband polarization at lob • This signature 
of Bloch oscillations (BO's) has been seen, not only in the 
THz radiation generated by the intraband polarization, 
but also in the time-resolved and time-integrated degen- 
erate four-wave mixing signals (TRFWM and TIFWM 
respectively). These features can all be qualitatively 
understood using perturbative approaches such as DCT 

S0. 

Recently, BO's have been observed in spectrally- 
resolved four- wave mixing (SRFWM) experiments |?3tllf|. 
where the peak energies of excitonic resonances were 
found to oscillate with the time delay, T21, at the Bloch 
frequency. Although it has been suggested that these os- 
cillations result from the influence of the excitonic intra- 
band dipole on the excitonic energies S 0, 0] , the pre- 
cise mechanism has never been clear and there has been 
some controversy regarding this interpretation [ll| . It 
is known that third-order DCT-type calculations do not 
yield significant peak oscillations [2j. In fact, as we shall 
show, peak oscillations similar to the experimentally- 
observed ones will only arise from a formalism that is 
infinite order in the optical field. 

In this letter, we present a new nonlinear response for- 
malism, valid to infinite order in the optical field, that 
combines the essential features of both the SBE's and 
DCT theory for asymmetric systems. Our approach re- 
lies on the fact that, in such systems, exciton-exciton 
interactions are dominated by a dipole-dipole interaction 
and phase space filling (PSF) effects are relatively in- 
significant for the range of carrier densities considered 
jl5j. In what follows, we present our theory and calcu- 
late the TIFWM and SRFWM signals for a BSSL. Wc 
obtain oscillations in the SRFWM peak energies that are 
in good agreement with experimental results HE3 . Fur- 
thermore, we find that the system is so nonlinear, that 
SRFWM and TIFWM results cannot be adequately de- 
scribed using a third order DCT approach, even at mod- 
erate densities (~3 x 10 9 /cm 2 ). 

We formulate our theory in the basis of the excitonic 
states of the BSSL in the presence of the applied dc field, 
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Edc- While recent work shows that treatment of PSF in 
an exciton basis is difficult and subtle 0, Q , these prob- 
lems do not arise for the densities considered here (lij . 
The excitons are characterized by the quantum numbers 
where /i describes the internal motion and K„ is 
the center of mass wave vector. This wave vector deter- 
mines the direction of the optical radiation emitted by the 
exciton. In degenerate FWM experiments, incident opti- 
cal pulses create excitons with wave vectors Ki = kiand 
K_i = \t2- These excitons produce an intrahand po- 
larization grating with wave vectors K 2 = ± (ki — k2) 
that scatter the optically-created excitons into wave vec- 
tors K3 = 2ki — k2 and K 3 = 2k2 — ki. Consider- 
ing all scattered excitons, the intraband grating is de- 
scribed by wave vectors K m = ?r (k x — k 2 ) (m even) 
and excitons are scattered into all wave vectors K„ = 
5 [(n + 1) ki — (n — 1) ka] {n odd) 0|- Henceforth, we 
refer to signals in the FWM direction, K_3, as FWM 
signals, although they are to infinite order in the optical 
field. 

Thus, for our BSSL, the Hamiltonian is given by Q 
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where V is the system volume, e is the permittivity of 
the BSSL, huj^ is the energy and B* n is the creation 
operator of the exciton (in the DC field) in the state 
n) (for n odd). The incident optical field is given by 
Re [E°p 1 (t)] = Re [e^ pt (t) e"^*] , where u c is the laser 
central frequency, £° pt (i) is the complex temporal en- 
velope, and (i) is an external THz field. The 
interband polarization with wave vector K n [n odd) is 

P l n "* er =^E M M /i 5 i,- n + h - c -i where M m is the exci " 
tonic interband dipole matrix element. Finally, the in- 
traband polarization with wave vector K m (m > and 
even) is P£'™ = $ E^, £„ G^B^ , n B^ <m+n , where 
pintra _ pintra* anc j sum over n is from — n to 
no — to by 2, and n Q is an odd, positive integer The 
index, n serves as a truncation point for the sum; if n a 
is infinite then the expression is exact. In actual calcu- 
lations we increase n Q until convergence is obtained. 

We use the Heisenberg equations to determine oper- 
ator dynamics and truncate the hierarchy of equations 
by factoring the three-exciton (six-particle) correlation 
functions into one- and two-exciton correlation functions. 
This factorization is very different than that used in the 
SBEs, as it retains the intra-excitonic electron-hole cor- 
relations and the long range exciton-exciton correlations. 
It has been shown to be very accurate for the calcula- 
tion of degenerate FWM signals to third order in BSSL's 
0. We account for dephasing and decoherence phe- 



nomenologically via the interband and intraband dephas- 
ing times, T M and T^ v respectively. The dynamical equa- 
tion for the (B^ n) thus becomes: 



dt 
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where E l ™4™ = -± (P™ tra ) + F^ z S m , is the total in- 
traband field with wave vector K m , and an equation sim- 
ilar to Eq. (J2J is used to obtain the {B* n B v< j) needed 
in the calculation of (P™ tra ). 

Note that in the above equations, the K = inter- 
nal intraband field interacts with the excitons in exactly 
the same way as an external THz field. The dc compo- 
nent of this field renormalizes the excitonic energy, fku^, 
while its THz component dynamically couples levels with 
different [l3|. As discussed below, it is largely the 
K = internal field that yields the peak oscillations in 
the SRFWM signal. A DCT calculation to any order 
will not yield these oscillations because the equation for 
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to order n) does not contain a term 
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proportional to E n nl \^>n-z 

We now present calculated results for the 
GaAs/Gap 7A.ln.3As superlattice used in recent ex- 
periments |9j, Il0| : the well widths and barrier widths 
are 67A and 17A respectively, the dc electric field is 
15 kV/cm, and there is no external THz field. The 
single-particle BO period of this system is tb = 328 fs, 
with a corresponding WSL energy spacing of fkus = 12.6 
meV. The system is excited by a Gaussian optical 
pulse with a temporal FWHM of 90 fs. The dephasing 



times are taken to be Ti, 



1.0 ps and 



1.5 ps, while the excitonic population 



decay time, TL,, is taken to be infinite ^(|. As in the 
experiments jjj llOj . the exciton areal density ranges 
from 10 9 to 10 10 /cm 2 . Only Is heavy-hole excitons are 
included, so that the single internal quantum number 
(i describes an exciton in which the average along-axis 
electron- hole separation is approximately fid Q. Thus, 
for the exciton, /1 plays the role that n does for the 
single-particle WSL states except that the excitonic 
WSL is not equally spaced Q. These states are 
calculated using the method described in Ref. 12]. The 
neglect of excited in-plane excitonic states has been 
shown to be justified for central laser frequencies below 
the u = WSL frequency [13 . 

Fig. 1 shows the TIFWM intensity as a function of 
delay time. The oscillations are due to BO's of the ex- 
citonic wavepackets 0, ^jj. The frequency of the os- 
cillations is not precisely given by wb for two reasons. 
First, the dc component of the K =0 internal THz field 
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rcnormalizes the excitonic energies by giving an effec- SRFWM polarization near lo = lo^ are approximately 
tive dc field of E^c + E ™ tra '' DC • because ^™ tra ' DC [ s 



antiparallel to ^dc f° r w c < w 0i the BO frequency is 
reduced, as observed. Second, excitonic effects alter the 
BO frequency by changing the WSL energy spacings as 
mentioned above. 

The large degree of nonlinearity in this system is seen 
in the TIFWM results obtained using different trunca- 
tion indices, n (dotted lines in Fig. 1). For no = 3, 
the TIFWM signal (K_3 direction) arises entirely from 
k2 excitons that scattered off the k2 — ki grating. As 
n is increased, the excitons can scatter out of K_3 into 
wave vectors K±„ for any odd n < hq and perhaps even 
scatter back into the FWM direction. Convergence is 
only obtained for n > 9 for this density, which demon- 
strates that scattering into wave vectors as large as K±g 
is clearly important and thus a third-order DCT ap- 
proach is not adequate. The multiple scattering and 
the high degree of nonlinearity is also evidenced in the 
fact that the intensities in the K_5 = 31t2 — 2ki and 
K_ 7 = 4k 2 — 3ki directions, also shown in Fig. 1, are 
comparable to the FWM intensity. 

In Fig. 2 we plot the SRFWM signal for a sequence of 
delay times. The spectral peaks are associated with dif- 
ferent excitonic states (as indicated in the figure) . There 
are three key features of these spectra that we draw 
attention to. First, the peaks do not occur precisely 
at the single-excition energies, w M ; this is due to the 
polarization-induced reduction of the applied dc field dis- 
cussed above. Second, some of the peaks deviate very 
considerably from a Lorentzian shape. Third, and of 
particular importance in this work, the peak positions 
depend significantly on the time delay, r 2 i. 

To see the peak oscillations more clearly, we plot in 
Fig. 3 the energy of the // = — 1 peak relative to the 
energy of the fj, = — 1 excition as a function of time delay 
for different exciton densities. As can be seen, and in 
agreement with experiment, the peaks' energies oscillate 
as a function of T21 with period given approximately by 
Tg. To understand qualitatively these oscillations, we 
consider a simplified model. We calculate the intraband 
fields to second order, use them in Eq. (J2J), and employ 
perturbation theory in all fields. The oscillations can 
then be largely understood as arising from quantum in- 
terference between multiple paths from the ground state 
to the excitonic state (/1, —3). In the inset to Fig. 2, 
we show the two lowest order paths. In the left-hand 
path, an exciton with wave vector K 1 = k 2 and energy 
fttOop is created by the second optical pulse; it then scat- 
ters off the K_2 grating into wave vector K_3, absorb- 
ing a photon of energy fkuf ntra ~ 0, ±Hujb- In the right 
hand path, the exciton additionally interacts with the 
spatially uniform field, E™ tra (t), absorbing a photon of 
energy huj® ntra ~ 0, ±Hujb before reaching the final state 
(p, —3). Including only the most important contribu- 
tions the left and right paths, their contributions to the 
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[cj - LOfj, + i/T mter ] 
(lo) E™ tra (lo b ) 

[lo - LO^ + i/Tinter] 



(3) 
(4) 



where E l Q ntra 



(u>b) is the Fourier component of El ntra (t) 
at the Bloch frequency. It is easily seen that interfer- 
ence between these two contributions yields peaks in the 
SRFWM spectra with energies that depend on the am- 
plitude and phase of E™ tra (ujb)- 

An approximate second order calculation for 
E l ntra (uj b ) yields 
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The oscillation of E l ntra 



(lob) with time delay is due to 
alternating constructive and destructive interference be- 
tween the Bloch oscillating intraband polarizations cre- 
ated by the first and second pulses. Thus, the periodic 
oscillations in the SRFWM peak positions arise due to 
the quantum interference between two-photon and three- 
photon processes. When the full equations are solved, 
one finds that higher-order nonlinearities play an impor- 
tant role, yielding an oscillation-amplitude that depends 
nonlinearly on exciton density. In the inset to Fig. 3, we 
show the amplitude of the peak oscillations (defined to be 
the energy difference between the first dip after t 2 \ = 
and the peak that follows) as a function of density for 
three different central laser frequencies. Note that at 
low densities, the oscillation amplitude is independent of 
density. In this regime, the peak oscillations are due to 
small effects that occur even to third order in the optical 
field. At higher densities the amplitudes increase su- 
perlinearly with density due to multiple scatterings from 
the intraband polarization grating. At the highest den- 
sity of 9.3 x 10 9 cm -2 , the sudden change in peak energy 
near T21 = 0.16 ps results from a splitting of the (J, = — 1 
peak, similar to that seen in Fig. 2 for the n = peak 
for r 2 i = 0.34 ps. 

Most of the features of the calculated peak oscillations 
agree with experiment: the phase, frequency, dependence 
on lo c , and la rge amplitude at T21 = 0, are all in general 
agreement [fl llOj . However, the amplitude of the peak 
oscillation is considerably smaller than that obtained ex- 
perimentally. For example, for lo c = loq — 2.27lob, and a 
density of 10 10 cm -2 , the calculated amplitude is roughly 
0.6 meV, while the experimental amplitude is approxi- 
mately 2.3 meV 0]. We believe that this difference is 
largely due to the effects of the screening of the exciton- 
exciton interactions via incoherent carriers. As discussed 
in Ref. [Io|. this plasma screening will considerably de- 
crease the dielectric constant at the Bloch frequency. We 
find, for example, that a decrease in e/e from 12.5 to 8 
can result in in increase in the the peak-oscillation am- 
plitude by a factor of 3 or more, bringing the calculated 
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and experimental results into quite good agreement. For 
simplicity, we have used the static value of 12.5, but plan 
to use a dynamic model in future calculations. 

In summary, we have developed a new infinite order ex- 
citonic system of equations for calculating the nonlinear 
optical response of asymmetric semiconductor nanostruc- 
tures. We have used these equations to reproduce for the 
first time the experimentally-observed [9j, LKJ oscillations 
in the SRFWM signals from BSSL's. 
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FIGURE CAPTIONS 

FIG. 1. The time-integrated intensity versus T21 for 
uj c = luq — 2.27lub, n — 13, and a density of 
6.36 x 10 9 cm 2 for the three different directions, 
K n indicated. For K_3, results for no = 3 (dash) 
and 5 (dot) are also plotted. 

FIG. 2. SRFWM intensity versus frequency at io c = 
ujq — 2.27u>b and a density of 9.3 x 10 9 cm 2 for a 
sequence of delay times. The inset is described in 
the text. 

FIG. 3. The fx = -1 SRFWM peak energy relative 
to fuv-i versus T21 at lj c = loq — 2.27lob for a series 
of densities. The curve offsets are a real effect due 
intraband renormalization of the DC field. The in- 
set shows peak oscillation amplitude as a function 
of density for (lu c — luq) /lub = —1.24 (squares), 
—2.27 (circles) and —2.83 (triangles). 
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